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The effects of parity and time reversal violating potential, in particular the tensor-pseudotensor 
electron nucleus interaction are studied. We establish that selected terms representing the interplay 
of these effects and the residual Coulomb interaction in the coupled-cluster method are equivalent 
to the coupled perturbed Hartree-Fock. We have shown that the normal CPHF diagrams have a 
one-one correspondance in the coupled-cluster theory, but the CPIfF pseudo diagrams are present in 
a subtle way. We have studied the pseudo diagrams in great detail and have shown explicitly their 
origin in coupled-cluster theory. This is demonstrated by considering the case of the permanent 
electric dipole moment of atomic Hg and our results are compared with the results of an earlier 
calculation. 
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I. INTRODUCTION 



The observation of a non-zero intrinsic electric dipole 
moment (EDM) of a non-degenerate quantum system is 
an evidence of parity(P) and time-reversal(T) symmetry 
violations [HH]. Among the two, T violation is of particu- 
lar interest as it is less understood and has important im- 
plications for physics beyond the standard model. In the 
experiments where CP violation has been observed so far 
[SlHlIS], T violation is inferred [6^ by invoking the CPT 
theorem. However, the observation of an EDM would 
be a direct evidence of T violation in nature. Atoms are 
suitable and promising candidates to measure permanent 
EDMs due to their sensitivity to the P and T violating 
phenomena in the nuclear (hadronic), electron-nucleus 
(semi-leptonic) and electron (leptonic) sectors [7]. In 
this paper, we study the atomic EDM arising from the 
tensor-pseudotensor electron-nucleus interactions, which 
is semileptonic in nature. The coupling constant Ct of 
this interaction is zero within the standard model but it 
is finite in some theories which are extensions of the stan- 
dard model [8^. The closed shell atoms are sensitive to 
this interaction due to its dependence on the nuclear spin 
and heavier atoms are more sensitive to this interaction 
since it scales as [5]. The closed shell atoms on which 
EDM experiments have been performed to date are ^^^Hg 
[inj and ^^^Xe [TTJ [T^l [13] and efforts are underway to 
improve the results. In addition, new experiments are 
planned for Yb[Tll [TB] and Ra fW . The ^^^Hg experi- 
mental data has been used to provide improved limits on 
important P,T odd coupling constants at the elementary 
particle level [TUl. These have been extracted by combin- 
ing atomic calculations [17] and experiments [10]. The 
coupled-perturbed Hartree Fock (CPHF) [18] effects are 
extremely important in the calculation of atomic proper- 



ties and it is particularly true in the case of the atomic 
EDMs [19]. It is important to analyze these effects in 
the framework of an all-order many-body method like 
the coupled-cluster method (20j, one of the most accurate 
many-body methods for the study of atomic properties 
[21] . The main thrust of this work is to demonstrate 
that all the CPHF effects are subsumed in the coupled- 
cluster method. This is an important step towards im- 
proving the accuracy of the existing atomic calculations 
of closed-shell atomic EDMs, which is necessary for ob- 
taining better limits on the P and T violating coupling 
constants. A similar study has been done for parity non- 
conservation effects in open shell atoms [22 . This pa- 
per is organized as follows : In Section |ll] we present 
the theoretical background, in Section [Hi] we discuss the 
coupled-cluster equations with/ without the T-PT inter- 
action, in Section IIVI we discuss the atomic EDM in the 



coupled-cluster and the CPHF framework, in Section |V] 
we discuss our results and finally we present our conclu- 
sions in Section IVTl 



II. THEORETICAL BACKGROUND 

A. Dirac-Fock equation 

For heavy atoms, the relativistic effects cannot be ne- 
glected and should be incorporated in the atomic Hamil- 
tonian. An approximate relativistic atomic Hamiltonian, 
appropriate for our calculations is the Dirac-Coulomb 
Hamiltonian i?DC- For ^ electron atomic system in 
atomic units 

N ^ I 
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where c is velocity of light, a and /3 are the Dirac ma- 
trices and fij is the separation between the ith and jth 
electrons. This Hamiltonian is an approximation to the 
general atomic Hamiltonian, which consists of additional 
terms arising from the interaction of electron spin with 
nuclear spin, interaction between spins of electrons, mag- 
netic interactions, etc. However, these interactions are 
neglected as their strength is negligible compared to the 
electron-electron Coulomb interaction and the nuclear 
potential energy. 

The single electron equations are obtained by approx- 
imating the two-electron Coulomb interaction by the 
Dirac- Fock central potential J/dfI^i)- Then 

N 

Hue = [ca^ ■ + Ac^ + VN{n) + UuFin)] + V^s, 

i 

(2) 

where the residual Coulomb interaction 

The residual Coulomb interaction embodies the non- 
central or correlation effects, which can be incorporated 
in the atomic theory calculations as a perturbation. The 
single electron wavefunctions satisfy the Schroedinger 
equation 

[ca .p + (3c' + F^(r) + UuF{r^] 1^-°) = e°lV'°>, (3) 

where j'i/'a) is the single electron wavefunction and a de- 
notes the quantum numbers which specify the wavefunc- 
tion uniquely and e° are the single electron energies. We 
can group the operators in the equation and rewrite the 
equation as 



(4) 



In the above equation t = ca. ■ p + (it? + V/v(r) 
and the Dirac-Fock potential 3°|^a) — 

E^T M\vW^m - {i^'tlvmm] , ^ere 
V — l/(ri — r2), ri being the position coordinate 
of the ith electron and the summation is over all the 
occupied orbitals in the reference state. If an atom has 
a non-zero EDM, then it is an indication of P and T 
violating interactions within the atom. In this paper, we 
consider the atomic EDM arising from the tensor-pseudo 
tensor P and T violating electron-nuclear interaction 



i 2V2CTGFl3a.- Ip{r), 



(5) 



where Gp and p{r) are the Fermi coupling constant and 
nuclear density respectively, (3 and a are the Dirac matri- 
ces and / is the nuclear spin. This interaction perturbes 
and the orbitals acquire an admixture from the op- 
posite parity orbitals. At the single electron level, these 
effects are incorporated in the CPHF calculations. 



B. Coupled-perturbed Hartree-Fock 

The introduction of the P and T violating interaction, 
/iT-PT, as a perturbation modifies the atomic Hamilto- 
nian. The corresponding single electron wavefunctions 
are the mixed parity states 



(6) 



where A is the perturbation parameter and is the 
first order correction, which is opposite in parity to ji/'a)- 
However, there is no first order energy correction as 
/iT-PT is parity odd. Then the perturbed Dirac-Fock 
equation is 



6=1 



'Y{%kb\v\'4>a)\i>b 



(7) 



6=1 



Selecting terms linear in A and rearranging we get the 
CPHF equation 

(/,0 + _ ^OJ 1^1^ ^ (-/.T-PT - g') Hi). (8) 

The perturbed Dirac-Fock operator 



b=l 

{^l\v\rl:l)m - {i^l\vm\^l)]. (9) 



In the present work, to solve Eq.([8|, we expand \'4'a) in 
terms of a complete set of unperturbed opposite parity 
orbitals. Then, \'tl>l) = Y.pGpa\'^l), where Cpa are the 
mixing coefficients. Substituting this in Eq.Q and pro- 
jecting by ('i/'pli we obtain a set of linear algebraic equa- 
tions 



Cpa ~ Eq) + E] l^^qabC'*^ + VpbaqCqh 
bq 



{p\h. 



T-PT 



la) = 0, 



(10) 



where Vpqab = {pq\v\ab) - {pq\v\ba) and Vpbaq = 
{pb\v\aq) — {pb\v\qa). Another approach to calculate 1-0^) 
is to solve Eq.(|8| self-consistently [T7]. Hereafter, for 
brevity, an orbital is represented as Eq.(10l can 
be written as the matrix equation 



AC 



(11) 



where Apa = Y^bqi^pqab + Vpbaq + (fip - Ea) SpqSab) and 
Bpa = (p| /iT-PT I a). This equation is solved iteratively 
starting with the zeroth order (in v) contribution 

^(0,1) ^ 
^pa 



gO _ gO 



(12) 
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(i) (ii) (iii) (iv) (V) 

FIG. 1; CPHF diagrams at zero and one order Coulomb interaction, v . The diagrams (ii,iii) are the pseudo CPHF diagrams 
and (iv,v) are the normal CPHF diagrams. The dotted Une is the Coulomb interaction and the line attached with Kl is the 
T-PT interaction. Pseudo diagrams are the diagrams with local energy denominators, note the direction of 6 and q orbitals. 



as the initial guess. The coefficients in the k*"^ iteration 
are 



Bpa J2bq 



VpqabC. 



(fc-1,1) 
qb 



VpbaqC. 



(fc-l,i: 
qb 



p 

(13) 

The superscripts in Cpa'^'' refer to the order of v and 
/iT-PT respectively. The diagrams arising from the above 
equation in zero and one order of v are shown in Figjl] 
In this paper, the diagrams of the second and the third 
terms in Eq.( 13 1 are referred to as the pseudo and normal 



diagrams respectively. 



III. COUPLED-CLUSTER EQUATIONS 



In the coupled cluster theory, the exact atomic state 



|*)-e^""|<i>o) 



(14) 

where T*^"^ is the cluster operator and |<i>o) is the ref- 
erence state. For the ground state, |<i>o) is the Slater- 
determinant of all the occupied orbitals. The cluster op- 
erator r'^' = J2i Ti^^ , where T^-^^ are the i-tuple exci- 
tation operators. The cluster amplitude equations are 
obtained from, after applying the operator e~"^' ' and 
projecting on excited states, the Schroedinger equation of 
l^-). Restricting to the approximation T^") = t[°^ +T^°\ 
the cluster operators are hence solutions of the equation 
((f>*|i?Ar|<i>o) = 0, where |<I>*) denotes singly and doubly 
excited states and respectively. For any oper- 

ator O, O = e~^'°' Oe^'°' is the dressed operator. It is to 
be noted that Hn is the normal ordered atomic Hamil- 
tonian, which is -ffoc in the present calculations. Let the 
-ffx-PT perturbed atomic state 



(15) 



where T^^^^ is the i/x-PT perturbed cluster operator and 
A as defined earlier, is the perturbation parameter. The 
perturbed coupled cluster equations are 



<i>*'| \Hn,T^'A 1*0) = -(<!>* 'I 



-T-^PT 
N 



(16) 



^T-PT 



where |<i>* ) are opposite in parity to |<i>o) and Hj^ 
is the normal ordered dressed _ffT-PT- Defining (O1O2} 



as the contraction and normal order of two operators Oi 
and O2, then 



[i/^,r(i)] = {HArr(i)}. 



(17) 



Consider T^^) = T^^\ the singles equation from Eq.p6| 
is 



(18) 



The cluster operator T^^^ — J2a p '^p'^ata *a 
sociated cluster amplitude. Retaining only hJ^^^^ from 

— T-PT 

H , in terms of matrix elements 



bq 



(19) 



Then the perturbed cluster amplitudes tP^^\ are solu- 
tions of the iterative equation 



JO) _ (0) ■ 



(20) 



This is equivalent to the Eq.(13l, the equation of the 



CPHF mixing coefficients without the pseudo diagrams. 
This formally establishes that the norma/ diagrams in the 
CPHF approach are equivalent to a subset of terms in the 
coupled-cluster theory. However, a similar comparison 
of the pseudo diagrams in the CPHF and the coupled- 
cluster theories is done later as it requires the dressed 
electric dipole operator D, defined in the next section 
through the atomic EDM. 



IV. ATOMIC EDM 
A. General expression 

The atomic EDM in the CPHF approximation is 
D,^J2 {a\d\p)ClT'^ + Cpa*^°"''\p\d\a), 

ap 

= 2 J2 c^pT''\ (21) 



ap 
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where d is the single particle electric dipole operator and 
the first superscript on the mixing coefficients refers to 
all order in v. The diagrams which contribute to Da are 
shown in Figj2] The mixing coefficients in Eq.(13l has 
V to all orders in the limit k —^ oo. Substituting the 
expression of CpT'^\ we have 



a\d\p) 



ap P 



bq 



-qb 



7r 

^pbaq'-^qb 



In coupled-cluster theory, the atomic EDM 



(22) 



(23) 



It should be noted that I?a has linked terms only the 
unlinked terms cancels the normalization factor. For the 



normal diagrams, the two terms in Eq.(23 1 are equivalent 
and arise from 13 (Ti^^). Here Tio^ is the contraction of 
T^'^^ with w-see Figlgiv) , (v) . 

B. Pseudo diagrams 

The diagrams in the Figj2] (ii) and (iii) are the sum of 
two many-body perturbation theory diagrams. We now 
show that, in coupled-cluster theory these are subsumed 

in Da = {(iyQllDTi'-^f^ +T^^^'' DT^"^]\^o)- Unlike in normal 



terms, here, Ti^g* = Ht-ptT^"' ■ The operator T^"' has 
correlation effects arising from the two-particle and two- 
hole as well as the other forms of w, to all orders. Hence, 
the CPHF pseudo diagrams do not have a one-one cor- 
respondence with the above CCEDM terms, rather the 
pseudo diagrams are part of terms in CCEDM. 

The algebraic expressions of the EDM diagrams shown 
in Fig. |2] (ii,iii) are. 



(0) 



and 



{a\D\p){pq\v\ab){b\hT^PT\q) 
(cp - fa) (e& - £«) 

{a\D\p){pq\v\ba){b\hT^PT\q) 

(Cp - ea) (Cb - Eg) 



(24) 



(25) 



respectively. The diagram in Fig |2|ii) is the sum of the 
two MBPT diagrams topologically equivalent to Fig. [3]- 
1(a) and 1(b), algebraically 



{a\D\p){pq\v\ab){b\hT^PT\q) 

+ <^b — — e<7 



(26) 



On simplification, this is same as Eq. ( 24 1 multiplied by a 
phase (-1). Similarly, the exchange diagram Fig. |2|iii) is 



the sum of the topologically equivalent MBPT diagrams 
of Fig j3] II (a) and (b), algebraically 



{a\D\p){pq\v\ba) /iedm | q) 
+ ^b — — e<7 



(27) 



which is equivalent to Eq. 25 apart from a phase — (- 
1 ) . To calculate the cluster amplitudes which has similar 

I 1 

correlation effects, retain i?T-PT and Ht-ptT^^ in Eq. 



( 18 1 . Then cluster amplitude equation is 



^ap 'l2bq ^pbaqta 



(28) 



where B'^p = {Ht~pt + HT--PTT!f'')ap- The cluster am- 
plitudes calculated from the above equation has effects 
of direct and exchange pseudo diagrams. 

To analyse the atomic EDM arising from the pseudo 
diagrams within the coupled-cluster theory, consider the 

two terms DTi''^ and t[^'^^ Dt!^°\ After Introducing 
complete set of eigen functions, the contribution to 
from these terms 



I 



(29) 



The first term on the right hand side, in terms of deter- 
minantal states, the contribution to the pseudo diagram 
is 



(30) 



this follows as D and t[^'^ are single particle operators. 

It is to be mentioned that, the {T^^)cs represents the 

I 1 

component of t[^^ arising from the term {HT-PTT^^)ap 



in Eq.p8|. Similarly, the second term in Eq.(29l is 



(<i>o|Ti('^Vs)«I^KD«|I^f I'&o). (31) 

ap,bq 



In this expresion, the pseudo diagram contributions are 
present in the component of arising from the term 
{HT-PT)ap in Eq.(28l. Using Slater-Condon rules, the 



matrix elements of single and two-particle operators be- 
tween determinantal states are: (<i>g|£)|<i>o) — {p\D\a) ; 

{Kl\4°^\'^o) = {pq\t2\ab)-{pq\t,\ba) and {<PI\D\<^>11) = 
{b\D\q). Then, 

Da = ^(p|i3|a)(a|(T(^))effb) + 

a,p 

hq^'^\a\D\p} [{pq\t2\ab) - {pq\t2\ba)]{32) 

ap^bq 
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(i) 





(iii) 





FIG. 2: CPHF diagrams contributing to EDM. (ii,iii) and (iv,v) are the EDM diagrams arising from pseudo and the normal 
CPHF diagrams respectively. The hnes with ® denote dipole operator. When Vcs is treated to all orders, the EDM interaction 
vertex is equivalent to the mixing coefficients in CPHF and singles cluster amplitude in CC. 
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1 cff 




+ 



(a) D(T<"W 



+ 




(b) T 




(I) 



(b) T . D T , 



(0) 



V. RESULTS 

A. Symmetrywise contribution 

From the expression of i7x-PTi the presence of nu- 
clear density p{r), Si/2 — P1/2 is expected to have the 
largest contribution. This is indeed observed in our cal- 
culations. The variation of Da for a small basis con- 
sisting of 68 Gaussian type orbitals [24j: (l-12)si/2, (2- 
13)pi/2,3/2, (3-10)^3/2,5/2, (4-7)/5/2,7/2 and (5-8)57/2,9/2, 
are given in Table |l] It lists the value of Da when or- 



FIG. 3: Diagrams contributing to EDM - Solid interaction 
lines in I(a)&(b), II(a)&(b) and III(a)&(b) represent the 
Coulomb interaction treated to all orders. The operator 

Ti^\g is a result of the contraction Tj^'r^'^' , which, when 
contracted with the induced dipole operator (D), gives the 
diagram contributing to Da- 



The the effects of the pseudo diagrams are distributed 
among various terms with the T2^'^ cluster amplitude. 
Hence, it is difficult to establish one-one correspondence 
between the CPHF pseudo diagrams and the correspond- 
ing diagrams in CCEDM in orders of v. This is a conse- 
quence of the structure of CCEDM and CPHF equations 
where the perturbed cluster amplitudes are computed us- 
ing the converged values of the Tg'^^ amplitudes, which 
treat the residual coulomb interaction to all orders. How- 
ever, the converged results which includes all orders of 
V and one order of _ffx-PT, where the sequence of the 
perturbations has all possible combinations, should be 
identical. It is possible to establish the equivalence by 
choosing only the two-particle two- hole terms of v in the 
2^2*^' equations. 

In this paper, we calculate the normal and the pseudo 
diagrams simultaneously, this couples the normal and 
pseudo contributions. This is evident from the Eq. 
( 13 1, where the CPHF coefficients are iterated with both 
the normal and the pseudo diagrams. But within the 
coupled-cluster theory, in particular CCEDM formalism, 
this inclusion is more subtle. This is due to the structure 
of Eq. ( 23 1 , where the converged T^^^ amplitudes contain 



the effects of the EDM from both the terms shown in Eq. 



TABLE I: Variation of Da with the inclusion of higher an- 
gular momentum virtual states. The first column implies that 
virtuals only upto the orbitals indicated have been included 
in the calculation, in addition to Si/2 and pi/2 symmetries. 
The normal and the two pseudo diagrams are calculated in- 
dependently. 

Virtual states EDM ( xIO"^^ e-m ) 

Normal (Normal+Pseudo) 



upto P3/2 
upto dr,/2 
upto fr/2 

upto g9/2 



-6.411 
-6.415 
-6.399 
-6.399 



-6.024 
-6.127 
-6.142 
-6.144 



bitals are added symmetrywise. According to the table 
Da , the contribution from the higher angular momentum 
d and / virtual orbitals are small and opposite in phase. 
Next, we consider an optimal basis set, with which we get 
converged Da, it consists 112 Gaussian type orbitals: (1- 



18)si/2, (2-18)p 



1/2,3/2, 



(3-16)4 



/2,5/2: 



(4-13)/; 



5/2,7/2 



and 



(5-10)37/2,9/2. The results from the optimal basis are 



TABLE IL Variation of Da with the inclusion of higher angu- 
lar momentum virtual states. The first column implies that 
virtuals only upto the orbitals indicated have been included in 
the calculation, in addition to Si/2 and pi/2 symmetries. The 
normal and the two pseudo diagrams have been calculated 
together. 

Virtual states EDM ( x 10"^^ e-m ) 



upto P3/2 
upto ds/2 
upto fr/2 

upto g9/2 



-5.83 
-5.90 
-6.75 
-6.75 



(281 



given in TablelTT] It lists Da arising from the normal and 



6 



pseudo, where these are calculated simultaneously. It can 
be seen that the contribution from the virtual orbitals 
enhances and changes the value —5.83 x lO^^^em to 
-6.75 X 10-22em. 

For a more detailed analysis, the dominant contribu- 
tions from 6si/2-Pi/2 ^^nd 6S1/2-P3/2 are listed in Table, 
nil The total contribution from the 6si/2~7^-Pi/2.3/2 from 



TABLE III: 



tI^'>DT^°'' (in units of 10" 



Dominant contributions to Da = 

~^^CTemaN) from the terms shown 
in Eq. [28] for np intermediate states calculated using the 
coupled-cluster theory for EDMs. 



Occ. 


np 




D 


Da 

T^'^D 


6Si/2 


6pi/2 


104.69 


0.872 


-1.010 


6Si/2 


7pi/2 


-254.88 


-1.821 


-5.139 


6Si/2 


8pi/2 


262.28 


1.388 


-4.032 


6Si/2 


9pi/2 


-202.37 


-0.344 


-0.771 


6Si/2 


lOpi/2 


-113.94 


0.068 


0.0858 


6Si/2 




-56.22 


0.692 


0.0431 


6Si/2 


6p3/2 


13.85 


0.995 


0.153 


6Si/2 


7p3/2 


-36.27 


-2.372 


0.953 


6Si/2 


8p3/2 


-36.80 


-2.211 


0.901 


6Si/2 


9p3/2 


15.91 


0.771 


0.0135 


Total 








-8.2026 



the term t[^^' DT.^ 
the term DT^'^^^ is 



n(0) 



is —.337 X 10 ^"^CTevacTM and from 
8.20 X 10-22 emCTCTAT. 



B. Cluster amplitudes and Da 

The calculated normal T^^"^ amplitudes are in excellent 
agreement with the corresponding CPHF mixing coeffi- 
cients. Calculating the normal and the pseudo diagrams 
together, the Da of atomic Hg is —6.75 x 10~^'^em. It is 
enhanced to —6.92 x 10~'^'^em when the two pseudo terms 
are calculated separately. A previous calculation [TT re- 
ported the CPHF Da of atomic Hg as —6.0 x 10~^^em. 
We attribute the difference of our result from the pre- 
vious calculation to the inclusion of correlation effects 
in CC, beyond those present in CPHF, which are dis- 
cussed in Section HVl The different numerical methods 
used can also contribute to the discrepancy, however this 
would be small. The different results, when the pseudo 
terms are calculated together and separately with the 
normal terms, is the effect of the coupling between the 
two terms. Comparing the contributions from the normal 
and the pseudo diagrams, the pseudo diagrams though 
important, are just ~ 9 % of the normal diagram con- 
tribution. The portion of the pseudo diagram contribu- 
tion is however dependent on the size of the basis set. 
For example, with the basis (l-14)si/2, (2-14)pi/2. 3/2,(3- 



12)^3/2,5/2, (4-8)/5/2,7/2 and (5-9)^7/2,9/2, the contribu- 
tion from the pseudo diagrams is 4 %. This indicates that 
the contribution from the pseudo diagrams increases till 
it converge. The phase of the normal diagrams is deter- 
mined by the phase of most domininant term, the Dirac- 
Fock contribution. For ^^^Hg this is negative. On the 
other hand, the phase of the pseudo diagrams cannot be 
ascertained easily. Unlike the normal terms, the lead- 
ing contribution from pseudo diagrams has T2^'^ which 



contributes to two terms 



Ht-ptT^°^ and DT: 



(0) 



The 



phases of the dominant Ht^ptT2^^ and DT^' contri- 
butions determines the overall phase of the pseudo di- 
agrams. That is, the relative phase of the normal and 
pseudo diagram is not a general trend. It depends on the 
phase of the of the dominant T^^ cluster amphtudes and 
hence it is atom specific. 



n(0) 



VI. CONCLUSION 

In this paper. We have numerically tested and demon- 
strated the inclusion of CPHF effects in coupled-cluster 
for atomic ^^^Hg. We have shown that there are cer- 
tain terms in the coupled-cluster theory for EDMs which 
are equivalent to the normal diagrams in the CPHF the- 
ory. This is demonstrated in the context of the EDM 
of i^^Hg arising from the T-PT electron-nucleus interac- 
tion, a property which is sensitive to the accuracy of the 
wavefunctions in the nuclear regions. The equivalence of 
the pseudo diagrams is more subtle and unlike normal 
cannot be shown explicitly. However, we identify the 
terms in coupled-cluster which correspond to the pseudo 
diagrams based on an analysis using many-body pertur- 
bation theory. The pseudo diagrams are the sum of two 
many-body perturbation theory diagrams [19 . Hence, in 
the coupled-cluster expression of Da in Eq.(23l, the di- 



rect and the conjugate terms when added give exactly the 
pseudo diagrams of CPHF. This shows that the coupled- 
cluster theory contains all the CPHF effects . The rela- 
tive phases of the pseudo and the normal diagrams are 
atom-specific and hence cannot be generalized. An opti- 
mal basis consisting of 112 orbitals gives converged result 
and is in good agreement with the result of Martensson- 
Pendrill [17] . We have studied and analysed in detail the 
various many-body effects that play an important role 
in EDM of atomic Hg and impact of coupling between 
normal and pseudo diagrams. 
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TABLE IV: Dominant contributions to Da{in units of 10 '^^CremaN) from the term DT2 



Occ. 


np 




(atomic units) 




^^CTemfTjv) 






Direct 


Exchange 


Direct 


Exchange 


6-51/2 


(yv i /■■} 


0.0926 


-1.077 


-0.011 


0.012 


Gsi/2 


7'D-\ r? 


-0.238 


0.269 


-0.067 


0.076 




8pi/2 


0.274 


-0.287 
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